THE THIRD MOMENT OF QUADRATIC DIRICHLET L-FUNCTIONS 

MATTHEW P. YOUNG 

Abstract. We study the third moment of quadratic Dirichlet i-functions, obtaining an 
error term of size 0{X^^^~^'^). 
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^ ■ 1. Introduction 

^ , Quadratic twists of L-functions have been extensively studied from many points of view. In 

this paper, we consider the problem of the third moment of quadratic Dirichlet L-functions. 
Soundararajan [S] was the first to find an asymptotic formula for this third moment, obtain- 
ing a power savings in the error term. Subsequently, Diaconu-Goldfeld-Hoffstein |DGHj used 
the multiple Dirichlet series method and obtained a superior error term. Our goal in this 
paper is to further refine the error term in this problem by adapting a technique introduced 
in 0. 
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Theorem 1.1. Let F be a smooth, compactly- supported function on the positive reals with 
support in a dyadic interval [X/2,3X], and satisfying 



^: (1.1) F(^)(x)«,X-^ J = 0,1,2,.... 

nr ! Then 



in 

en 

O ■ where the star indicates the sum is over squarefree integers, and P{x) is a certain degree 6 

polynomial. 



:i.2) 5^*L(l/2,X8dfi^(rf) = $^*P(logrf)F(rf) + 0(X3/^+^ 

(d,2)=l (d,2)=l 



For the smoothed sum as above, Soundararajan [S] previously obtained an error of size 
^ I X^/*"*"^ (for this, see the final displayed equation on p. 487 and optimally choose Y = X^^^). 

Using the multiple Dirichlet series method, [DGH] obtained an error of size X^/^"*"^ for the 
smoothed mean value. 

We shall present the full details of the proof of Theorem 1 1 . H but the technique can certainly 
be generalized to cover other cases. For instance, one can find an asymptotic formula for the 
second moment of the same family, perhaps with an error term of size 0{X^^'^'^^). It would 
be interesting to study the first moment of quadratic twists of a GL^ automorphic form such 
as the symmetric-square lift of a holomorphic modular form on the full modular group (we 
mention this case since the Ramanujan conjecture is known for such forms which avoids some 
potential pitfalls). This case should be largely similar to the third moment in this paper, but 
there is an extra difficulty since one cannot exploit the factorization of the L-functions and 
use Heath-Brown's quadratic large sieve |H-B] (see comments following Theorem 13 . 1 1 below) . 



This material is based upon work supported by the National Science Foundation under agreement No. 
DMS-0758235. Any opinions, findings and conclusions or recommendations expressed in this material are 
those of the authors and do not necessarily reflect the views of the National Science Foundation. 
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We stress that the condition that d is squarefree is a substantial difficulty (improperly seen 
as a simple technicality). For instance, with F as in Theorem II. ![ compare the difficulty in 
the following two sums: 

(1.3) A = Y,nd\ B = Y,*nd). 

Poisson summation quickly shows A = F{0) + 0{X~'-'') for any large C > 0. On the other 
hand, the standard Perron-type formula approach quickly shows B = jjJ-+0{X^^'^), and any 

improvement on the exponent 1/2 would give a quasi- Riemann hypothesis, that is, a zero- 
free region for the Riemann zeta function in Re(s) > 1 — S for some 6 > 0. The underlying 
reason for this is that the generating function for the squarefree numbers is ({s)/((2s). 

The method of proof of Theorem 11.11 uses the main idea of ^ which is a kind of recursive 
argument that efficiently treats the squarefree condition. In effect, there is almost no cost 
to handling squarefree numbers compared to all integers (up to a barrier at improving on 
an error of size 0{X^^'^^^); our method unfortunately does not lead to a quasi- Riemann 
hypothesis!). 

Here we briefly sketch the method. We start with an approximate functional equation for 
the central value (actually we consider values slightly shifted from the central point). The 
basic idea is to use Mobius inversion on the sum over d to remove the squarefree condition, 
and then to use Poisson summation. The Mobius procedure effectively shortens the sum 

over d: say we write 2, fi^) = '^a^^('^) '^df('^'^d), then for large values of a the sum 
over d would be quite short indeed, and then Poisson summation might not be a wise choice 
as the dual sum would be longer than the original sum. The idea is to treat small and 
large a differently. For a small one uses Poisson summation; for a large one could use the 
trivial bound as in [5]. However, for a large one can obtain a more precise result by writing 
d — )■ b^d where the new d is squarefree. Then one arrives at something similar to where 
one started, and rather than bounding this expression trivially one can plug in an already- 
obtained asymptotic formula. Doing so leads to certain partial main terms that turn out 
to combine with other main terms from a small. This was the main new idea appearing 
in 0. The idea is simple yet carrying out the details involves proving elaborate (perhaps 
miraculous?) combinatorial type identities for the various main terms that arise in different 
ways. 

We attempted to perform the same technique on the third moment in this paper, and 
the method largely succeeds though there are some new difficulties. One difficulty is that 
the main terms do not combine quite as nicely as in 0; apparently we are missing some 
harmonics. This seems to be related to the problem of getting an asymptotic formula for 
the fourth moment of this family (where if one uses Poisson summation and collects all the 
main terms arising in similar ways as in this paper, then one does not obtain all the main 
terms predicted by [CFKRS] ) . To get around this issue of missing terms (which are not main 
terms but rather unwieldly expressions that should combine and simplify), we noticed that 
a particular choice of a weight function in the approximate functional equation can cause 
some of these terms to vanish. 

Obtaining an error term better than 0(X^/^) is of great interest for a few reasons. For 
one, it is apparently related to the subconvexity problem for a GL^ L-function twisted by 
a quadratic character. The basic idea is that if the central values were non-negative, and 
the weight function in Theorem 11.11 were allowed to be supported in a short interval (as in 
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[S]), then we could drop all but one term in the moment and obtain an upper bound of 
size 0{X^^'^~^) for that L- value, which would be a subconvexity bound as the conductor 
is X X^ (note that recently Blomer [B] obtained subconvexity for quadratic twists of a 
symmetric-square lift of a Maass or holomorphic form of full level, though not using the 
above approach). Although we cannot rigorously claim that improving the error term in this 
way implies subconvexity, it seems that the issues are related. 

The combinatorial identities in this paper are quite elaborate. This is a general phenome- 
non observed in various examples yet not well understood. In some previous works we found 
some techniques that simplify the proofs of these identities, which involve showing that two 
very different arithmetical Dirichlet series turn out to be the same, at least for special values 
of the parameters. Because of the sizes of the expressions involved here, we resorted to 
computer verification of certain identities. We stress that the remarkable simplifications in 
this paper give very strong evidence to their correctness; our experience shows that even a 
tiny error completely destroys the simplifications. 

As in any proof by induction, it is extremely beneficial to have the inductive hypothesis 
provided ahead of time. This is accomplished by the general moment conjectures of |CFKRS] ; 
their method leads to the form of the main term very quickly and with minimal effort. 

It is interesting to compare how our approach is related to the multiple Dirichlet series 
method used in |DGHj . To this end, we briefly summarize their method in simplified terms. 
One might wish to study the two-variable Dirichlet series 

Z{s,w)= ^ . 

d>l, odd 

By the Perron formula method, the meromorphic continuation of this function in terms of w, 
with s = 1/2, is intimately related to the error term in (11. 2p . In fact. Theorem 11.11 gives the 
meromorphic continuation of Z(l/2, w) to Re(w) > 3/4 with a singularity a.t w = 1 only (we 
include a brief sketch of this fact in the next paragraph). One approach to studying Z{s, w), 
initially in its region of absolute convergence, would proceed by writing the Dirichlet series 
expansion for L{s, XsdY, using Mobius inversion to remove the condition that d is squarefree, 
and then applying Poisson summation to the sum over d. Heuristically, this should lead to a 
connection between Z{s, w) and Z{s, 1 — w), but unfortunately the Poisson formula does not 
lead to exactly the same Dirichlet series (we see this in our work in Section 15.11 and Lemma 
15.21 below). Instead, the method of |DGHj is to construct a different auxiliary Dirichlet 
series say Z'{s, w) that does satisfy nice functional equations but differs from Z{s, w) in that 
it involves non-fundamental discriminants and corresponding correction factors at "bad" 
primes in the Euler products. The functional equations of Z' make it a pleasant function to 
understand, and with a kind of sieving argument this information can be transferred to Z 
(and hence to the moment on the left hand side of (II. 2p ). 

To see that Z{l/2,w) has the desired meromorphic continuation, let F be a function 
satisfying the conditions in Theorem 11.11 with X = 1, and such that F{x) > for all x 
(but not identically zero). With F denoting the Mellin transform of F, we have F{w)d~^ = 
J^ x~'^F{d/x)^. Note that the real part of F{w) is always positive (this is direct from the 
definition) and hence this function has no zeros. Then 



Fiw)Zil/2,w)= / x-'"( Y,* Lil/2,xsd)'Fid/x)) 



d>l, odd 
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Applying (11. 2p . we obtain a main term plus an error term. The main term can be explicitly 
computed which gives the meromorphic continuation of this term, while the error term has 
absolute convergence for Re(w) > 3/4 and hence leads to the desired analytic continuation. 
It would be desirable to extend the class of allowable weight functions in Theorem 11.11 in 
order to obtain non- vanishing results in short intervals, or to un-smooth the sum in (11. 2p . 
To obtain strong results in this direction, it would be desirable to study moments analogous 
to (11.21) but at points s with potentially large imaginary part. The multiple Dirichlet series 
approach is well-suited to keeping track of this s-dependence. As an aside, we mention that 
un-smoothing (11.21) is somewhat subtle because the central values are not known to be non- 
negative; as a result, the error term as stated in Theorem 2 of [S] is not justified without 
this assumption. It would also be valuable to understand the other classes of fundamental 
discriminants in place of those of the form 8d with d odd squarefree and positive; in principle, 
the different cases should be similar to each other. 

1.1. Acknowledgements. I thank Brian Conrey, David Farmer, Jeff Hoffstein, Mike Ru- 
binstein, and K. Soundararajan for discussions. 

2. Tools 

In this section we quote some results we use throughout the paper. For d a fundamental 
discriminant, let Xd be the corresponding primitive quadratic Dirichlet character. We shall 
work with discriminants of the form 8d where d is odd, squarefree, and positive, so that 
Xsdin) = (^) for n odd is an even, primitive character of conductor Sd. 

2.1. Functional equation. The functional equation of the Dirichlet L- function associated 
to xsd is 

— J r(0L(s,X8d) = A(l-s,X8d). 
In its asymmetric form it reads 



(2.2) Hs,xu) = X(s)Hl-s,Xad), X(l/2 + u ' ^'' 



" r {^) 



IT J r I V2+" 



2.2. Approximate functional equation. Suppose that a,/3,7, called the "shift param- 
eters", are small complex numbers. We shall suppose initially that each A G {a,/3,7} 
lies in a punctured rectangle of the form |Re(A)| < Ci/logX, |Im(A)| < C2X^ minus 
|R'e(A)| < ci/(21ogX), |Im(A)| < (c2/2)X^ where the q (depending on the choice of A) 
are chosen so that the domain corresponding to each shift parameter has distance ^ 1/X^ 
from the other two shift parameters. Later we can relax this condition. 

Proposition 2.1 (Approximate functional equation). Let G{s) he an entire, even function 
with rapid decay in the strip \Re{s)\ < 10. Then for Xsd os above, 

Xsdin)aa,i3,^in) / n 



(2.3) L(| + a, XsdW^ + /3, XsdWl + 7, Xsd) = J2 ^''.'^''^'"' v^,,,, ( 

n n2 



c/3/2 






n 



77,2 
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where notation is as follows: 

(2-4) K,/3,7(a;) = TT- / -^9a,/s,yis)x-'ds; 

zm j(^i^ s 



^'■'^ ^w^) = i-j TTiT^TTwf 77F^^ ' 



2 



ro,/3,7 = rQ,r/3r^, where 



(2.6) 



r 



1 

2~" 



^/ -r. / k+a^ 



r 



2 



(2-7) (^a,/3,7H = Yl 



a-'^b-^c- 



abc=n 



Remark 2.2. We shall choose G to vanish at the poles o/C(l + 2a + 2s), C(l + 2/3 + 2s), 
etc., and to be divisible by a large collection of values of ^{s) = s(l — s)7r~*/^r(s/2)C(s), 
including ^(4 + 2a + 2(3 + 4s), ^(2 + 2a + 27 + 4s), etc. (here we have listed the first couple 
terms in the numerator and denominator, respectively, of (14.321) with a,/3,7, replaced by 
a + s, /3 + s,7 + sj. More precisely, these zeta factors occur in the analysis of the arithmetical 
factor Aa+s,i3+s,-y+s dcfincd by (I2.13P and developed much further in Lemma \4.1\ below. For 



a given e > small we choose G{s) to vanish at the poles of all the C 's which occur in 



Lemma 4-1 as numerators (i.e., with da.b,c > ^), and also to be divisible by all the C's which 



occur in Lemma 4-1 as denominators (i.e., with da.b,c < 0/, when obtaining the meromorphic 
continuation of Aa+s,i3+s,"f+s to Re{s) > — | + e. 

We furthermore suppose by a symmetrization argument that G{s) is symmetric under 
any permutation o/{a,/3,7}, and under switching any a,(3,'y with its negative, and under 
switching s with —s. Then by scaling we can ensure G{0) = 1; Assuming the shift parameters 
lie in the punctured rectangles as above then in terms of the shift parameters we have G{s) ^ 
X'. 

Making such a choice for G remarkably simpUfies certain later computations. The point is 
that we wish to use a contour shift argument to analyze certain integrals in the development 
of the moment under consideration. In so doing one would cross many poles arising from 
the arithmetical factors. However, the contribution from these residues would involve values 
of G{s) at points other than s = and since G can be chosen from a wide class of functions, 
it is apparently unlikely that these residues can persist as main terms. By having G vanish 
at these points we see a priori that these terms do not survive in the final answer. In |Y] . 
we did not move the contours past such poles but instead matched up integrals to form the 
simplifications. One might hope to see a similar matching here but in fact it seems some 
of the terms are missing, that is, the method does not capture all these terms. One would 
naturally speculate that these missing terms arise from an incomplete analysis of the new 
mean value arising after Poisson summation (i.e., in (I5.35P below). 
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2.3. Poisson summation. We now quote Soundararajan's result 

Lemma 2.3 ([S]). Let F be a smooth function with compact support on IR+, and suppose 
that n is an odd integer. Then 



(d,2)=l ^ ^ ^ ^ fcez 



where 



'-) '^'<"'-(¥-(v)^) E ©<^ 

^ ^ '' ^ a (mod n) ^ 

and 

/oo 
(cos(27rxy) + sin(27rxi/))F(x)rfx. 
-oo 

The Gauss-type sum is calculated exactly with the following 

Lemma 2.4 ([S])- //"i andn are relatively prime oddintegers, thenGkimn) = Gk{m)Gk{n), 
and if p" is the largest power of p dividing k (setting a = oo if k = 0), then 



(2.11) Gkip" 



0, if (3 < a is odd, 

(pip^), if f3 < a is even, 

—p"', if P = a + 1 is even, 

(^)p"v^, z//3 = a + l tsodd, 

0, z//3>a + 2. 



2.4. The |CFKRS] conjecture. Suppose F{d) is as in Theorem ll.il Then 

Conjecture 2.5 ( |CFKRS] ). Suppose a,/3,'j lie in the rectangle \Re{s)\ < j^^, 1/772(5)1 < 
X^ . Then if I = I1I2 is odd with li squarefree and I2 a square, then 

(2.12) Y,* ^(i + «' Xm)L{\ + (3, xm)L{\ + 7, Xsd)x^d{l)F{d) 



(d,2)=l 

where 

(2.13) A.,,,(o= $: ^-Miii!^ n(i+^'')-' 

(n,2)=l p|ra/i'2 

= (2(1 + 2a)C2(l + 2/3)C2(l + 27)C2(1 + « + /3)C2(1 + « + 7)C2(1 + /? + 7)5a,/3,-y(/), 

where -Bq,./3,^ /ios an absolutely convergent Euler product for the parameters in a neighborhood 
of the origin. Furthermore, the meaning of TJl^^' ^ is simply T^^V oVt^ , where 5i = if 
ei = +1, and 6i = 1 if ei = —1. 
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The conjecture is derived from the following orthogonality relation 



i\-i 



(2.14) Yl* Xs,{m)F{d) ~ ^ 1[{1 + p' 

(d,2)=l ^^^ ' p\m 

for m an odd square, and is o{X) otherwise (for fixed m). This relation appears in a slightly 
different form in |CFKRS] but can be directly calculated using Poisson summation along the 
lines of the calculations in Section 15.11 

The expected value of / is somewhat controversial. The five authors |CFKRS] conjectured 
that / = 1/2 is allowable, while |DGHj predict the existence of a main term with / = 3/4; the 
subtlety here is that |DGHj analyze the moments involving non- fundamental discriminants 
and correction factors at bad primes. A sieving argument is required to convert between the 
different moments. Q. Zhang [Z] provided further evidence for / = 3/4 and in fact gave a 
prediction for the numerical value of the constant in the lower-order term of size X^'^. The 
constant is ~ —.2 making the lower-order term difficult to detect with numerics; however, M. 
Alderson and M. Rubinstein |ARj recently performed extensive calculations which perhaps 
show modest agreement with this (numerically small) lower-order term. 

3. Outline of the method 
Our recursive approach to the problem takes the form 
Theorem 3.1. // Conjecture \2.5\ is true with f > 3/4, then it is true for f replaced by 

3 _|_ J"! 

4 "^ 2/ • 

The case / = 1 is "trivial" in the sense that the error term is larger than the main 
term, but quite nontrivial in that it uses the quadratic large sieve of Heath-Brown |H-Bj . 
Actually, we could start with / = 7/4 which is an immediate consequence of the convexity 
bound (however, we emphasize that we still require the quadratic large sieve in the proof 
of Theorem 13.11 though perhaps with extra work it could be removed). In any event, the 
sequence /i = 1, /2 = 7/8, /s = 23/28, etc. leads to / = 3/4 as allowable (one can check 
that the sequence of /„'s is decreasing and bounded below by 3/4 and hence has a limit 
which is easily checked to be 3/4). 

3.1. Dissection. Before embarking on the details of the proof, we shall give an overview of 
the whole argument, deferring proofs to later sections. First we need some notation. Let 
M{1) = Ma^p^^{l) be the moment on the left hand side of (12.120 . Using the "two-piece" 
approximate functional equation from Proposition 12.11 write 

(3.1) M„,;3,^(/) = Mi(/) + M_i(/), 

respective to the two sums in the approximate functional equation. That is, 

(3.2) M, = y: m E ^"^"^^r'^^'V ..,, ^ ^ 



rf3/2 



(d,2)=l n n2 

and 

(3.3) AU = r„,,, J:' d-'-^--F(d) J: 2«£NK^.t2M v_„,_,,_, ( " 

(d,2)=l n n2 
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The expressions are similar enough that we may focus upon Mi and then easily deduce 
analogous formulas for M_i. Indeed, M_i is the same as Mi after swapping a and —a, /3 
and — /3, 7 and —7, replacing F{x) by F_q,^_^^_^(x) = x~°'~^~'^F{x), and multiplying by 
Fq,^^^, in that order. 

As explained in the introduction, the plan is to use Mobius inversion on the sum over d 
and to treat the resulting sum in two ways. We begin by removing the condition that d is 
squarefree, getting 

(3.4) Ml = 2^ fi{a) 2^ F{da ) }_^ j ^.,^3,7 ( ^^ 1 . 

(a,20=l (rf,2)=l {n,2a)=l n2 \l y' / 

Now we separate the terms with a <Y and with a > Y [Y a. parameter to be chosen later), 
writing Mi = Mn + M^, respectively. We similarly write M_i = M_7v + M^r. 

3.2. Overview of Mj^j. To evaluate Mjv, we use Poisson summation. Lemma 12.31 on the 
sum over rf, getting a sum over fc, say. The term A; = gives a certain main term, which 
we denote M^ik = 0). As Soundararajan realized [5], the terms fc 7^ should be separated 
into squares and non-squares. The squares give three more main terms which we denote 
M^{k = □,«), Mj^ik = n,/?), M]s[{k = 0,7). Evaluating these "off-diagonal" main terms 
is much more subtle than for k = and is accomplished in Section [5] below. The work with 
Mtv is summarized with the following 

Lemma 3.2. For the special choice of G{s) given by Remark \2.^ we have 
(3.5) 
Mat = M^(A; = 0) + M^(fc = D, a) + M^(A; = D, /3) + M^ik = □, 7) + 0{YI^'^+'X^I^+'). 

The A; 7^ terms can be naturally expressed as a certain contour integral; moving contours 
to the left picks up poles giving the three "off-diagonal" main terms. Bounding the contour 
integral on the new lines of integration gives the error term. 

3.3. Overview of Mr. The evaluation of Mr is in some sense simpler than for M^ as it 
is purely combinatorial. However, the expressions become somewhat complicated; here we 
give a sketchy argument that we make rigorous in Section H] below. Recall 

(3.6) Mr= ^ fi{a) ^ F{da ) ^ j \4,^,^ ( 1 . 

(a.20=l (d,2)=l (n.2a)=l W^ \\ I / 

a>Y 

Now apply the change of variables d — )■ b'^d with the new d squarefree, to get 

(3.7) M« = $: ,(a) ^ $:• F(d(am E ""'""'?'^''''"V „... (77^^) ■ 

{a,2l)=l {b,2l)=l{d,2)=l (n,2ab)=l 77,2 \\\ ) ) / 

a>Y 

Using the definition of V as an integral representation (12. 4p . we get that the inner sum over 
n above is 

in,2ab)=l n2 '' (^) 

This sum over n almost produces a product of L-functions but the coprimality restriction 
(n, 2ab) = 1 slightly perturbs it. Ignoring this annoyance for the present discussion, we see 
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that Mfi should be related to 

(3-9) XI ^'^"^ XI ^ / ^— ^",/3,7(*)^c.+s,/3+.,7+s(0c^S, 

(a,2/)=l {b,2/)=l "^^^ 

a>Y 

where M^^s^^j^s^^^siV) is the same moment with which we started but with a different weight 
function having support with integers of size X/{ab)'^. Then we can apply Conjecture 12.51 to 
this sum, expressing it as a sum of eight main terms (one for each choice of ei, e2, £3 G {±1}), 
plus an error term. 

After making these arguments rigorous (treating the coprimality restriction {n,2ab) = 1, 
etc.), these ideas lead to 

Lemma 3.3. // Coniecture \2.5\ holds with a parameter f > 1/2, then 
(3.10) Mn= Yl MR(ei,e2,e3) + 0(V/^^), 

where M/j(ei, €2, 63) is defined below as 04. 8p . 



3.4. How the main terms combine. Of the eight main terms appearing in Lemma 13.31 
only four turn out to contribute significantly. We suspect that finding an error better than 
X^^^ would (at least) require dealing with all eight of these terms. The smaller ones are 
bounded with the following 

Lemma 3.4. // exactly two of the e^ 's are —1, then for a special choice of G{s) described in 
Remark \2.2\ we have 

(3.11) Mniei, €2, €3) < YX^/\lXy, 
and furthermore, 

(3.12) Mij(-l,-l,-l)<X3/^(/X)^ 

This dispenses with four of the terms on the right hand side of (I3.10p . leaving four that 
do contribute. The remaining ones combine in a very pleasant way with the four main 
terms of Lemma [3.21 The term Mi^{k = 0) combines with M^{1, 1, 1), while Mi^{k = D, *) 
with * = a, /3, 7 combines with M/j(— 1, 1, 1), Mji{l, —1, 1), Mji{l, 1, —1), respectively. The 
combination is as pleasant as one could hope; each corresponding pair of M^ and M^j 
give almost identical expressions, the only difference being Mn has a < Y while M^j has 
a > y , so that combining them simply removes this restriction on the sum over a altogether. 
This matching of M^ and M^i is the only truly difficult part of this work, as it requires 
substantial calculation to see that these very different expressions actually agree. Having 
done the simpler case of the first moment |Yj, it was easier to predict that terms should 
combine in this way. 

Each of the four combined pairs M^v and Mji is given as a certain contour integral. By 
shifting the contour to the left, one crosses a pole at s = only; the residue at this point 
gives one of the eight main terms on the right hand side of (I2.12p . The integral along the 
new contour is estimated with absolute values. The end result is 

Lemma 3.5. For a special choice of G{s) described in Remark WTB, we have 

(3.13) M^(fc = 0) + M^,(l, 1, 1) = A^^p^,{l) f^^^ + 0(X^/^+T), 

2(,2l2jVH 
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and furthermore, 

(3.14) M,,ik = D, a) + Mni-l, 1, 1) = A„,,/3,,(/)r^ 5/) ~ '^ + ©(X^/^+^f 



and similarly for the other two terms (with a replaced by either /3 or 'y). 

By combining Lemmas 13. 2^ 13. 3^ 13.41 and 13. 5^ we get that Mi equals the sum of the four 
main terms on the right hand side of fl2.12p . namely those with at most one of the e^ = —1, 
plus an error of size 

(3.15) « /l/2+.x3/4+.y + lU2+e^^_ 

We then obtain an asymptotic for M„i by the procedure described following (13.31) . An 
examination of the form of the main terms on the right hand side of (I2.12p shows that M_i 
accounts for the remaining four main terms in (I2.12p , namely those with at least two of the 
Si = —1, plus the same error as given in (13.150 . Thus, taken together, these results show 

Lemma 3.6. // Coniecture \2.5\ holds with some f > 1/2, then 

(3.16) Ml + M_i = M.T. + 0(/^/2+"X3/'+"y) + 0{-^j^), 
where M.T. is the main term on the right hand side of (I2.12p . 

Choosing Y = X '^f gives Theorem 13. 1[ It has been so far left implicit that the error 
terms in the above lemmas are uniform with respect to the shift parameters lying in the 
appropriate punctured rectangles described in Section 12.21 In Lemma 13.61 we can claim the 
same error term for shift parameters in the non-punctured rectangles |Re(A)| < Ci/logX, 
|Im(A) < C2X^ by the following reasoning. We note that Mi + M_i is analytic in terms of 
the shift parameters in this larger domain; by work of |CFKRS] . the main term in (I2.12p 
has this same property. Therefore, so does the error term. The maximum modulus principle 
shows the error term is maximized in the punctured rectangle domain where we initially 
proved the result, and so the uniformity is extended to the larger domain as claimed. 

We prove Lemma [3.31 in Section WA[ Lemma [3.41 in Section l42l Lemma [3.21 in Section [5^ 
and Lemma 13.51 in Section Taking these four results for granted, we already showed that 
Lemma [XS] is valid, and hence Theorem 13. II holds. The five authors |CFKRS] clearly explain 
how to let the shift parameters tend to zero and hence deduce Theorem 11.11 from Theorem 

EH 

4. Calculating with Mr 
In this section we obtain Lemma 13.31 and Lemma 13.41 
4.1. Proof of Lemma 13.31 
Proof. Our first subsidiary goal is to find the rigorous analog of (13. 9p . This formula is 



,4.1) M„^ Y. "(«) E E 4^1^1^5^/, E>«'<'^'^)-('--3) 

(a,2Z)=l (6,2/)=lri,r2,r3lab rf T^ r^ "^ (^)(d,2)=l 

a>Y 

G(s) 
L(i + a + s, X8d)L{l + /3 + s, X8d)^(| + 7 + «, X8d){da^bY'/^{rir2r3)-'^-^go,^f3^^{s)ds. 
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To prove this, we pick up with (13. 7p . recalhng (13. 8 P is the sum over n. Then note 

{n,2ab)=l p\2ab F F F 

= L{a + s,X8d)L{P + s,X8d)L{-f + s,X8d) 2^ ^+, ff+, ^+, • 

, T-\ To Try 

ri,r2,r3\ab i z o 

Inserting (14. 2 p into (13. 8p . and then into (13. 7p . gives (I4.ip . as desired. We were able to move 
the hne of integration to a = e since the Dirichlet L-functions have no poles. We take 
e = l/logX. 

Now we proceed directly to the proof of Lemma 13.31 Note that the inner sum over d 
appearing in (14. ip . namely 

(4.3) J2* X8d{lrir2r3)L{l + a + s,X8d)Lil + (3 + s,X8d)L{l + 'y + s,X8d){da%y^^F{da%^), 

{d,2)=l 

is of the form Ma+s,i3+s,'y+s{irir2r3) , but with a new weight function with smaller support 
{d-X/{abf). Write 

(4.4) {da%y/^F{da%^) = F,_^.,,^,,{d), 

where Fi,.y{x) = {xyYF{xy). 

Next we apply Conjecture 12.51 to this inner sum over d. Technically, we need to first 
truncate the s-integral so that |Im(s)| < (log(X/a^6^))^. When (a6)^ < X^~^, then the ex- 
ponential decay of the integrand shows that the error incurred by this truncation is negligible 
(<^ X~^^^ , say). Otherwise, when (a6)^ > X^~^ then the sum over d is almost bounded so 
that the convexity bound L{l/2 + a + s, xsd) ^ ((1 + kDMI)"^^"^ is sufficient to show these 
terms give 0(X2+^). After plugging in Conjecture 12.51 to the truncated integral, then the 
same argument works to extend the integral back to the whole vertical line, without intro- 
ducing a new error. By this procedure, we get Mpt as the sum of eight main terms plus an 
error of size 

(4.5) «$^|Ma)lE E / Ti 'V ^rir^ra --^ « v^TTiTn^'- 

a>Y ib,2)=lr,,r,,r,\ah (^^2^3)^ V^^/ ^^ 

This is the error term appearing in Lemma 13.31 

As for the main terms, we have by a direct application of Conjecture 12.51 that 

(Aa\ ij ( \ ^^ ( \ ^^ ^^ /^(n)/^(r2)/i(r3) 1 1 

(4.6) MR(ei,e2,e3) = 2^ ^l[a) ^ ^ 



{a,2l)=l 
a>Y 



^.^. . .^+° i+/3^i+7 ^(/nrarg)* 2(2(2) 

{b,2l)=lri,r2,r3\ab '1 '2 '3 V V J^ ^ •^y ^ \ / 



— I v4ei(a+s),e2(/3+s),e3(7+s)(^n?^2?"3)rj-|j'^}3+^_^+^ 



^3./2;a2fe2(l - 6i{a + s) - 52(/3 + s) - ^3(7 + s))- ga,l5,^,{s)ds. 

[rir2r3Y s 



12 MATTHEW P. YOUNG 

Here we use the notation m* = mi, where m = 77111712 where mi is squarefree and 7712 is a 
square. We can simphfy the Melhn transform a bit by noting 

~ f°° -A rlr ~ 

(4.7) F3,/2;,2,2(m) = / {xa'}?)'^F{xa'\?)x^— = (a6)"2«F(f + u). 

Jo X ' 

As shorthand, let w = 1 — 6i{a + s) — S2{f3 + s) — (^3(7 + s). Then 

(4.8) Mniei,e2,e,) = -^ ^1 ^'^''^^J ^itAs,,^s^9a,,,,is)nf + ^) 



(a,2l) 
a>Y 



El Y^ /^(^l)/^('^2)/^(^3) 1 . n \^ 

(b,2i)=l ^ ' ri,r2,r3|afe 'i '2 '3 V V J^ ^ •^/ 



These are the main terms appearing in Lemma 13. 3[ completing the proof. D 



4.2. Proof of Lemma 13.41 Here we prove Lemma [3.41 We shall need an auxiliary result 
that we first motivate. By a symmetry argument, for the proof of (13. lip it suffices to consider 
the case ei = e2 = — 1, £3 = 1. In this case, we have 

(4.9) M^(-l,-l,l) = -i- Yl M«)^/ T^+sTp+s^ga,p,,is)Fa-a-^-l) 

^^^ ' (a,2l)=l ^^ "^(^) ^ 



{a,2l) 
a>Y 



v^ 1 v^ /^(n)/^(r2)/i(r3) 1 ^ r, X, 

(6,2Z)=1 ^ ' ri,r2,r3|a6 'l '2 '3 VVi^J/ 

The important feature here is that the weight function satisfies F{1 — a — /3 — |) <^ X ~ 2+^^ 
where a = Re(s) so that moving the contour to the right gives a saving. The hitch in 
carrying this out is getting an analytic continuation of the integrand for sufficiently large a. 
For (T large the sum over a will not converge absolutely; to get around this, we extend a to 
all positive integers, and subtract the contribution from a < Y, writing Mr(— 1, —1, 1) = 
M'(-l, -1, 1) - M"(-l, -1, 1), respectively. 

First we simplify M'(— 1, —1, 1); by grouping ab into a new variable, say c, then the sum 
over a becomes Ylaicf^i^)^ whence c = 1 is the only term that does not vanish. Then b = 1, 
and Ti = r2 = r^ = 1, so that 

(4.10) M'(-l, -1, 1) = -^^ / r^^sTp^s^9a,p,,i^)F{l - « - /3 - f ) 

1 

-A-a-s,-fi-s,^+s\^-)d'S. 
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As for M"(— 1, —1, 1), we sum over b to get 

(4.11) M"(-l,-l,l) = -i- J2 '"^«)^/ ^c.+sT^+s^9a,fs,,is)F{l-a-P-l) 



(a,2l) 
a<Y 



1 >^ /i(ri)/i(r2)/i(r3) f {a,[ri,r2,rs]) 



2-2a-2/3-4s 



^2-2^-2/3-4. ^ !+„+, i+;3+, 1+^+, ^ [ri,r2,r3 

(rir2r3,2«)=l ' 1 '2 '3 

C2;(2 - 2a - 2/3 - 4s) ^ ^=A_a_^ _/3_a,^+a(/rir2r3)rfg. 

V{lrir2r3)* 

We shall estimate both (14.101) and (14. lip in essentially the same way, namely by moving 
the contour of integration to a = ^ — e and bounding everything with absolute values. We 
shall show that the sums over the Tj converge absolutely. In this way we reduce the problem 
to obtaining the analytic behavior of the arithmetical factor to this region. This behavior is 
detailed with the following 

Lemma 4.1. Let I = I1I2 and Aa,i3,j{l) be as in Conjecture \2.5[ Then Aa,i3,j{l) has a 
meromorphic continuation to Re{a) , Re{f3) , Re{'y) > —\. Precisely, for any positive integer 
M there exist integers da^b,c (possibly negative or zero) such that 



(4.12) v4«,/3,7(0 = C'a,/3,^(/) Yl ({a + b + c + 2aa + 2bP + 2c-fy- 



b,c 



a+b+c<AI~l 
a,b,c>0 



where for any 6 > 0, Ca^^'^il) is given by an absolutely convergent Euler product in the region 
Re{a) , Re{(3) , Re{'j) > ^^^ + ^- Furthermore, in this region Ca./3,'y{l) satisfies the bound 



(4.13) a,/3,7(/) < Vhl'. 

Remark. If a + 6 + c = 1 then da,b,c = 1? which justifies the representation on the second 
line of (12131) . 

Before proving Lemma [4.11 we conclude the proof of Lemma \3M We require the estimate 

(4.14) F(l-a-/3-f)<X'-f, 

and that G{s) is prescribed according to Remark 12.21 Thus by moving the contour to 
0" = I — 5, we have 

(4.15) M'(-i, -1, 1) < /i"^/f ^'^x^-^/^+^ 

and 
(4.16) 



M"(-l,-l,l)«X3/^-^^/^ya^^ y (^^i^2r3)-|Mri)Mr2)Mr3)| Aa,[ri,r2,r3])\ 

^ ^tl (rir2r3)i-^ V [n,r2,r3] / 



4<5 



a<Y ri,r2,rs 



Here e > is arbitrarily small. Using the trivial inequality (a, [ri,r2,r3\) < a and noting 
that the sums over ri, r2, r^ converge absolutely for any 6 > 0, (provided e < 6/2, say). Thus 
by taking 6 > arbitrarily small, we get 

(4.17) M"(-l, -1, 1) < YX'^/\lXy. 



14 MATTHEW P. YOUNG 

These bounds furnish the proof of (13. lip . The case of fl3.12p is similar but easier. In this 
case, 



(4.18) 






G{s) 



Mb(-1,-1,-1) = — — - > /i(a)— - / r„+,,^+,,^+,^-^^„,^,^(s)F(l-a-/3-7-f) 



(a,20= 



>r^ 1 >r^ /i(ri)/i(r2)/i(r3) A^^-s -/3-.,-7^s(^^i^2^3) , 

^ (a6)2(i-"-/3-7-3s) 2^ i+„+, 1+/3+. 1+7+s yri;^;^^^^ 

(b,2l)=l ri,r2,rs\ab 'l '2 '3 \/ \ i z ^j 

In this case we simply move Re(s) = a to a = | — 5 for 5 small. The sums over a and b 
converge absolutely so an easy argument gives (I3!l2p . 



Proof of Lemma 4-i Recall 

(4.19) A^,,,,ii) = j2 ""^f^"'^ n(i + p^y'^ 

(n,2)=l p\nl 

which is initially defined for Re(a, (3, 7) > 0. The function A has an Euler product represen- 
tation, say A = YlpAp- Write /i = IlpP'''' say. Then 

(4.20) ^.=i: "-''-y' n(i+^-)-'. 

3=0 P q\lpi 

The (1 + Q'"^)"^ term is mildly annoying and is benign with respect to proving Lemma [4. H 
so we first make a simple approximation to remove this factor. For p\ 21, Ap has the form 






/ 2j'\ 

(4.21) A, = 1 + (1 + p^'r' J2 """'f = 1 + (1 + P~r\-i + 4 

say, where 

00 
(4.22) A'p = J2 



(ya,fia{P^^) 



j=0 ^ 

Note that for Re(a), Re(/3), Re(7) > -i, 

(4.23) A'p = l + 0{p~\p-^'' + p-^^ + p-2^ + p-"-'^ + p-°-^ + p-''-^)). 
Thus 

(4.24) ^ = 1 + lizi±^ = 1 + o(p"2(p-2" +p"2^ +p-2^ +p"""^ +p-"-^ +p"^-^)), 
Ap p zip 

so that Y[p\2i ~^ is absolutely convergent in any domain of the form Re(a), Re(/9), Re(7) > 
— I + 5 with 5 > fixed. Thus W^^x Ap has the meromorphy properties required by Lemma 
14.11 if and only if Y[p\2i ^'p does. (We shall treat the primes dividing / later). 
Note that 

(4.25) cr„,^,^(/)= Yl P"'^"-"^-^^. 

a+b+c=k 
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Thus for p \ 21, 

j=0 ^ a,b,c>0 P 






Let X = p 2 '^^ y = p 2 /^^ and z = p 2 t. Then with this notation, 

(4.27) a; = ^ ^"'^'7 ^ = E ^(1 + (-1)'^^''"'^)^^^, 

j=0 ^ a,fe,c>0 

which simphfies as 

1 1 1 1 + xy + xz + yz 



(4.28) Ap 2^(l-x)(l-2/)(l-z) ' {l+x){l + y){l + z)' (1 - x2)(l - y2)(i _ ^2) • 
Let -D(x, y, z) = {\ — x^)(l — 2/^)(l — z^){l — xy){\ — xz){l — yz), whence 

where 

(4.30) N{x, y,z) = 1 — x^y"^ — x^z^ — y^z^ — x^yz — xy^z — xyz^ + (degree 6 and higher). 
Furthermore, N[x, y, z) G Z,[x, y, z]. Then we can write 

(4.31) Ar(x, y, z) = {\- x^y'){\ - x^z^){l - y^z^){l - x^yz){l - xy^z){l - xyz^)N2{x, y, z) 

where N2{x, y,z) = 1 + (degree 6 and higher). It is clear that this process can be continued 
indefinitely to obtain A' as a ratio of products of polynomials of the form (1 — x'^y^z^) with 
i+j + k less than any given bound, say B, times a polynomial of the form 1 plus a symmetric 
polynomial of degree > B + 2. In terms of Y[p\2i ^p' ^^^^ gives a factorization in terms of 
ratios of the zeta function, of the form 

(21(1 + 2a)C2i{l + 2/3X21(1 + 2^X21 



(4.32) n4 



^^^ C2i{2 + 2a + 2^X21(2 + 2a + 2^X21(2 + 2/3 + 27) 

C2/(l + « + f3X2i{l + « + 7)C2/(1 + /3 + 7) 



(21(2 + 2« + /3 + ^X2i{2 + « + 2/3 + ^X2i{2 + a + /3 + 27) ■ ■ ■ ' 

where the dots indicate the contribution from the terms corresponding to the degree 6 
factors of N2{x,y, z). The degree 6 factors give an absolutely convergent Euler product for 
Re(a), Re(/3), Re(7) > —^. By continuing this procedure, we continue extracting ratios of 
zeta functions, each of the form C2i{M + aa + b(3 + 07) with a + b + c = 2M; by taking such 
terms corresponding to degree < M — 1 polynomials of x, y, z, we get absolute convergence 
for the "remainder" term provided Re(a), Re(/3), Re(7) > ^^^ which as M — > 00 leads to 
Re(Q!), Re(/3), Re(7) > — |, as desired. 

As for p\l, first notice that extending the product in (I4.32p to p\l gives the ratio of zeta 
functions as in (I4.12p . and that this finite product is holomorphic and bounded by /^ in the 
region Re(a), Re(/3), Re(7) > — | + 5. 

To complete the proof of Lemma |4?T| we need to examine the behavior of YlpU ^p- ■^'^^ Pl^ 

(4.33) A, = {i+p~r'J2^^ — ^• 

j=o ^ 
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Here Ip = or 1. The work above easily handles the case /p = 0, so suppose Ip = 1. Then 

(4.34) J2^I^lM}l—^= J2 p~-'--bf5-c-y-^i±^±^^ 

j=0 ^ a+b+c=l (mod 2) 

which with the notation x, y, z as before, equals 
(4.35) 

V/^Z^2^ V ; ; y ^■^2Vl-x 1-y 1 



a,b,c 



(1 - x)(l - i/)(l - ^) (1 + x){l + y){l + zy 



Simplifying this, we get that f l4.34p equals 

(A'ia\ x + y + z + xyz 

^^■^^^ ^ (l-x^)(l-y^)(l-.^) - 

For Re(a,/9,7) > — | + 5 > — ^ this is analytic and satisfies the bound 

(4.37) < |p-°| + \p-^\ + \p~^\ < p'^^ 

with an implied constant depending on 6. Thus J^ .^ Ap is holomorphic in the desired region, 
and satisfies 

(4.38) J]Ap«/Jr, 

p\i 

as desired. This completes the proof of Lemma 14.11 D 

5. Development of Mn 

5.1. Poisson summation. We shall use Poisson summation to analyze M^. We perform 
these calculations here. Recall that Mn is given by (13. 4p but with the additional truncation 
condition a <Y imposed. We write 

(5.1) M.. E M».) E (^)^ E (^) ^C'"^)'^- (i^ 

(a,2Z)=l (n,2a)=l ^ ^ n2 (d,2)=l ^ ^ \V ^ 



By Poisson summation. Lemma [2.31 this is 
(5,2) 



(a,2Z)=l (n,2a)=l ^ ^ n2 fcgZ 



(a,2Z)=l (n,2a) 

a<y 



where 



(5.3) H{y)= {cos + sm){27ixy)F{xa^)Va,i3,^ i , ^ .3^^ j dx. 

We find an alternate Mellin-type expression for H with the following. 
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Lemma 5.1. Suppose that Cg^Cu are real numbers such that ^ — c„ > and Cg > 0. Then 
for y ^ 0, we have 



1 / 1 V /• f ~. .^(S^ ^ -2 



(^•^) ^(^) - ^ l« j i, 1, ^(^ - ")^«-.(') i^ 1 " 



2 



) S TT I S 

r(Y -M)(cos + sgn(i/)sm)(-(y - u))dsdu. 



Furthermore, for c^ > 0, we /iawe 

1 1 /■ G{s)^ ,„^_.£;,, , 35 

Proof. First suppose 2/ 7^ 0. We begin by writing 



(5.5) H{Q) = -— / ^^g^^p^^[s)n-^F{l + -)ds. 

a^ ZTTZ ./r^ -1 s 2 



(5.6) F(xa2) = — f F{l + u){xa^)-^-''du, 

valid for any c„ G M, and the Mellin formula for Va^p^^y, i.e., (12. 4p . Then we change variables 
X — >■ a;/(27r|?/|) to obtain 

(5.7) iJ(y) = l r(cos + sgn(y)sin)(x)f-i.') / / F(l + «)^^.,^,,(.) 



a^ \ ^ " 33 „ , , da; 



n " — — - X 2 ^dsdu — . 

\2'K\y\J X 

Next we reverse the order of integrals in order to perform the x-integral first. We temporarily 
make the assumption < ^ — c„ < 1 for convergence. We argue as in my previous paper 
to justify the reversal of integrals. Then using 

f dx TTin 

(5.8) / CS{x)x'"— = Tiw)CS{—), 

Jo X 2 

where CS stands for either cos or sin, and < Re(w) < 1, we obtain (15. 4p . 

The case |/ = is much easier. We simply change variables x — )■ x/a^, apply the definition 
(12. 4p . reverse the orders of integration, and evaluate the x-integral. D 



Let M^ik = 0) denote the contribution to Mn from k = 0. Applying Lemma [5. II to (15. 2p . 
we have 
(5.9) 

M„(i = 0) = l Y. ^^[ ^S.,,{s)hl + 'i) E '-^£f^^ds. 

(a,2/)=l "^^^^l {n,2a)=l ^2 

a<Y 



By Lemma [2. 4[ we have Go{nl) = (j){nl) if nl is a square, and is zero otherwise. This means 
we can write n — )■ /in^, and then we obtain 

(5.10) M^{k = Q) = \ Y. ^^ [ ^9a,(s,-y{s)Fil + ^)D^{k = 0, s)ds, 



(a,2l)=l " '" -^i^) 

a<Y 
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where 

(5.11) D,{k = 0-s)= V ^-^-M0[M_ 

Now consider the terms with k j^ 0, which we denote M^^k ^ 0). Suppose e G {±} is the 
sign of k. Then we obtain Mjqik 7^ 0) = M^(k 7^ 0) + M^(A; 7^ 0), where by applying Lemma 
IS.ll to (15.21) and reversing the orders of summation and integration (supposing 0^ = 0^^ 3, 
say), we have 

,5.., MMM0)4 E ^(^)7ot''^^"'^^"-<^'(v *" 

a<Y 

Next we argue with inclusion-exclusion to treat the term (—1)'^. Suppose that f{k) is some 
function such that /(4/c) = A~^f{k)] in our case, we have f{k) = Gek{nl)/\k\^~'^, and 
z = Y — M. Note that nl is odd and hence G^^ekip-l) = Gek{nl). Then we write k = kik^ with 
ki squarefree and separate ki odd and even. Then we have 



(5.13) Y.^-iff{k) = E* E /(^i^2) + E* E (-i)'v(fcifc 

fc=l fci even A:2 fci odd k2=l 



(5-14) = E* E /(^i^') + E*(2 E f^^^^^l) - E /(^i^'))- 

fci even fc2=l fci odd fc2=l ^2=1 

Then using f{4:k) = 4~^f{k), this becomes 

oo oo oo 

(5.15) $^(-1) V(A;) = (2^-^^ - 1) E* E /(^i^2) + E* E /(^i^^2) 

fe=l fei odd fc2=l fci even fc2=l 

Applying the decomposition (I5.15P to M^{k 7^ 0), we obtain 

(5.16) M^(MO) = i Y. ^(E*-^i(^'^'^i'^)+ E*-^2(5,n,A;i,0), 



2 ^^ a^ 

{a,2l)=l ki odd fci even 

a<Y 



where 



3s 



(5.17) M,{s,u,k,J)={±-X [ [ F{l + u)^g^,,,,{s)(^) (2^-^^^'^ - i) 



r(— — M)(cos+esin)(-(— — u))Jekii'is -2u,\ + u- ^)dsdu, 



and where 

(5-18) J...(.,-)= E E^^4^ ^^' 

(n,2a)=lfc2>l ^ 
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We suppress the dependence on /, a, a, f3, and 7 in our notation for J. The formula for 
M2{s,u,kiJ) is identical to f l5.17p except the factor (2^"^'^"'"^" — 1) is omitted. 

5.2. Continuation of J. Here we develop the analytic properties of Jekiiv,w). 

Lemma 5.2. Suppose that e = ±1, {1, 2a) = 1, ki is squarefree, and J^kA'^iW) is given 
initially by flS.lSp for say Re{v) > 2 and Re{w) > 2. Then for any 6 > 0, Jeki{v,w) has a 
meromorphic continuation to Re{w) > 6 and Re{v) > 2, provided a,^,^ are small enough 
compared to 6. Furthermore, in this region we have 

(5.19) Jefci(f , W) = L2al{\ + W + a, Xek^)L2al{\ + W + (3 , Xek^)L2al{\ + If + 7, Xekjhkj_{v, w) , 

where Xeki (p) = ( ~ ) ? ^^^ subscript on L2ai means the Euler factors dividing 2al are re- 
moved, and where Ieki{v,w) is analytic in this domain and satisfies the bound 

(5.20) hkA^,w)<^s,el-^^'. 

Proof. We begin by observing that J has an Euler product as we now explain. Lemma 12.41 
immediately shows the multiplicativity in terms of n. It follows from the definition of Gk{n) 
(and a change of variables) that if k is multiplied by a square coprime to n then the value 
of Gk{n) is unchanged. 

According to the Euler product, we write J^ki{v,w) = Y[p'J<iki;p{'v,w). In the easiest case 
p\2a, we have 



(5-21) J,krA^,w) = J2^, = {^ 

If p I 2a, and if say p^p \\l, then 



p'"r\ 

j=0 ^ 



n=0 fc2=0 ^ ^ 

Suppose p \ ki. In view of Lemma 12.41 there are two classes of terms where the Gauss 
sum does not vanish, one of which is for n + Ip even and n + Ip < 2k2, and the other is for 
n + Ip = 2/^2 + 1. Thus 

r523) / (vw)= V V ^^A.(p") 0(p"^'-) , (f) V ^"./^.^(P''^^"'^) 

^ '' ^eA.-i;plty, ^ J /_^ ^ ^nw+kiV jfi+lp "^ /^ Z^ ^{2k2+l~lp)w+k2V ' 

n=lp {-mod 2) 2k2>n+ly ^ ^ V^ 2fc2>/p-l 

Executing the sum over k2 in the first sum above, we obtain 
(5.24) 

-JekyA'^^) {^ P ) 2^ +(Il+il)„ pn+l, ^ /p 2^ {2k2 + l-l,)w+k2V ' 

n=lp (mod 2) P ^ 2k2>lp — l 

li p\ki then the calculation is similar except we use the third line of (12. lip . The two classes 
of solutions come from n + Ip even and n + Ip < 2^2 + 1, and from n + Ip = 2k2 + 2. The 
former class gives an identical contribution to the previous case since for n + Ip even the 
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condition n + Ip < 2/^2 + 1 is equivalent to n + Ip < 2k2. Therefore in the case p\ki we have 
(5.25) 

n=«p (mod 2) P ^ ^ 2fc2+2>ip ^ 

A;2>0 

Now we analyze the above representations in closer detail. Suppose p\ 2akil. Then 
(5.26) 

(— ) 

Here e is small and accounts for the shift parameters. If Re(w) > 2 and Re(t(;) > 6, then 
(5.27) JekrA^,w) = [1 + (£^)p-i-(p-" +p-/^ +p-7) + 0{p-'+% 

which we write as 

The point is that the error term, when multiplied over all primes p\ 2alki is 05(1). 
Now consider the case p\ 2a/, p\ki. A short calculation shows 

(5.29) JekrA^^w) = (1 + Oip-'-^^+') + 0(p-2+-)), 

which pleasantly agrees with (I5.28P in case p\ki. 

Now suppose that p f 2aki, p\l, so /p G {1, 2}. If /p = 1 then 

(5.30) Jek,A^,w) = Oip"""-^^') + 0(p-^) 
Similarly, if /p = 2 then 

(5.31) Jek,A^,w) = Oip~n + 0(p"5"-"-+^). 

For each value of Ip, then, we have Jekup{v,w) ^ p-^v/"^ jn Re(t>) > 2, Re(ii') > 5. 
Finally, consider the case p f 2a, p|/ci, p\l. If /p = 1 then 

(5.32) J,fc,;p(t;, u;) = ©(p-'^-^'^^) + 0(j9-i-"'+^), 
while if /p = 2 then 

(5.33) J,fc,;p(i;, I/;) = 0(j>"^) + 0[p-^). 

For both values of /p, we have J^k^piv, w) ^ p-'-p/'^ in Re(f ) > 2, Re(M;) > 5. 

Gathering all these results completes the proof. D 

5.3. Proof of Lemma 13.21 We now have the ingredients in place to prove Lemma [321 For 
notational simplicity, consider the contribution from ki odd in 05. 161) . namely 

(5.34) I Y: ^E*-^^(^'^'^1'0- 

{a,2l)=l fci odd 

a<Y 
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The contribution from /c2 even follows similar lines and we omit its discussion. We write 

(5.35) Mi{s,u,h,l)=(^] I J H{s,u)(^y J,,^(3s-2M,i + n-f) 

s s s 

L2al{l + U- - + a, Xeki)L2al{^ + U - - + /3 , Xefci)^2a/(1 + M - - + 7, Xekjdsdu, 



where 



,<^(^)„ r.\^-^+urol-3s+2u iN-n/3s ^^^^_ , ..^N^^T^Ss 



(5.36) H{s,u) = F{l + u)^g^^^^,{s)7r--+-{2'-'^+'--l)T{--u){cos + sm){-{-~u)). 

Note that H has rapid decay in s and u in any vertical strips. We initially have Cg = c^ = 3, 
say. Next we move the contours (simultaneously) to c^ = c^ = | + e, in order to not cross 
any poles and to remain in a region where J^ki is known to be analytic. Next we move c„ to 
— ^ + e; in doing so we cross poles of the Dirichlet L-functions at m = | — a, w = | — /3, and 
M = I — 7 for e = fci = 1 only. We denote the contribution to M^ from these three residues 
as M]si{k = D,a), M^ik = D,/?), and M]^{k = 0,7), respectively. We shall develop these 
terms further in Section 16.21 and now focus on the error term coming from the new contour. 
The quadratic large sieve inequality of Heath-Brown |H-Bj shows 



(5.37) 5^* \L2al{(r + ^t)\'<^{alyK'+%l + \t\) 



l+e 

5 



K<ki<2K 



for I < o" < 1. Thus the sum over ki squarefree converges absolutely along these lines of 
integration. Furthermore, notice that with these choices of Cu and Cg that 

(5.38) His,u) « X3/4+^(l + |s|)-i°"'(l + \u\y'\ 
Hence the contribution to Mjv from these error terms is 

(5.39) < Yl ci-\la^Y+n-'-+'X^/^+' < /i/2+.^^3/4+.^ 

a<Y 

This is precisely what we needed to show for Lemma 13.21 It seems plausible that one 
could replace the application of Heath-Brown's quadratic large sieve with some of the work 
developed here since to a first approximation we obtain a third moment of quadratic Dirichlet 
L-functions in the above analysis. Strictly speaking, this is problematic because of the "bad" 
Euler factors and the presence of the factor I^^^ which does not match our intial setup. 
For future reference, we record the following expression for M^lk = D, a): 

(5.40) 

{a,2l)=l ^'^'''> ^ ^ 



a<Y 



r(s + a)(cos + sin)( — (s + a))Res^^i_„Ji(2s + 2oL,vj)ds. 



6. Proof of Lemma EIS] 
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6.1. Part 1. Here we prove that fl3.13p holds. We use the expression (IS.lOp . and compare 
it to 

(6.1) M^(l,l,l) = l Yl ^^ f ^a,a-^)Fa + ^)DR{l,l,l;s)ds, 

a>Y 

where 

(ao\ n n 1 1 ^ ^ \^ ^ \^ /i(ri)/x (r2)/i(r3) 



(■ (2) ^-^ h^ ^-^ l/2+a+s 1/2+/3+S 1/2+7+s 

^^ ' (b,2/)=l ri,r2,r3|a6 'l '2 '3 



(n,2)=l V 1 2 J/ / p|nirir2r3 



Using a computer algebra package, we shall check that in fact 

(6.3) D^(A; = 0;s) = D«(l,l,l;s). 

Notice that each Dirichlet series above has an Euler product and hence it suffices to check 
each Euler factor. For p\ al, we have that the Euler factor at p for D^ik = 0; s) is 



-aa-bp-c-y 
j>l ^ a+b+c=2j 



4-/3- 



-x'^y'z", 



Set X = p 2 '^ ^^ y = p 2 p ^ ^ z = p 2 T'', and define 

(6.5) T(x,y,^)= ^ ^^^^ x-^i/V, t/(x,|/,^)= 5^ ^ 

a,b,c>0 a,6,c>0 

which evaluate as rational functions in x, ?/, 2;. Then (16.41) takes the form 

(6.6) l + (l-p-i)(-l + T(x,y,2)). 

On the other hand, we compute the Euler factor at p for -D_r(1, 1, 1; s) ioi p\ al as 

(6.7) l + {l+p-^)-\-l + Q{x,y,z)), 

where Q corresponds to the same sum but with the annoying factor (1 + p~^)~^ removed, 
and a direct calculation gives 

(6.8) Q{x,y,z) = {l~p-^)[T + -—^ {r _ ^^^y ^ z)U + {xy + xz + yz)T - xyzU)]. 

p^(l — p~ ) 

Now we obtain an explicit rational function representation for (16.71) and a computer quickly 
verifies that these are identical. Similar computations cover the cases p\a and p\l (the final 
case p = 2 is trivial, both sides being 1), giving (16. 3p . 

Thus we obtain that 
(6.9) 

M;v(A; = 0) + M^(l,l,l) = l V ^-^ / ^^g^,p,,{s)F{l + ^-^)D^{k = ^-s)ds. 
Noting that 

{a,2ln)=l ^^ ' p\nl 
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and that 



(6.11) E TTw n(i + p~r' = ^A^,s,,,s,,,sii), 



(n,2)=l V^l"" ^- p\nl 



(n,2)=l 

we have 
(6.12) 

M^{k = 0) + M«(l, 1, 1) = -_i_=^ y ^^„^^_^(s)F(l + ^)A^+s,p+s,,+s{l)ds. 

By Lemma SrU and Remark l2.2l we can move the contour of integration to — 1/2 + e crossing 
a pole at s = only in the process. The residue at s = is as desired, and the error term is 
of size 0(X^/^"'"^/^), also as desired. 

6.2. Part 2. Here we prove (13.141) : this is much more intricate than proving (I3.13p . We 
begin with an Archimedean-type identity. 

Lemma 6.1. Let u be a complex number. Then 

(6.13) {2'-'- - l)(cos + sin)(^«)7r-"r(«) = 21.^^^^, 

where recall r„ is defined by (12. 6p . 

Proof. We first use the functional equation for the Riemann zeta function in the form 

(6.14) 7r-"r(M)C(2M) = 7r-^+"r(| - n)C(l - 2u), 
giving that the left hand side of (I6.13P is 

(6.15) (2i-2"-l)(cos + sin)(^«)7r-i+«r(i-t.)^^i^. 

By examining the Euler product, we observe that 

(6.16) (2i-2« _ 1)^(1 _ 2u) = 2^-2«^2(i _ 2u), 

so that the left hand side of (I6.13P becomes 

a-2«/„„^ , ^,„^,7^ 1 (2(1 -2m) 



(6.17) 2^"^"(cos + sin)(-«)7r-^+"r(i-«) ^^^^^ 

Next we use the trigonometric identity 

(6.18) cos(^) + sin(^) = V2cos(j-^), 
to get that the LHS of fl6A3D is 

(6.19) 2f-- cos(^(l - n))vr-i+"r(l - ^)^^|^. 
Next we use the gamma function identity 

(6.20) 7r-h^~''cosi^v)T{v^- ^^2) 



2^'-'^' r(if^)' 
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with V = ^ — u, getting now that fl6.13p is 



i-3..«r(V)C2(i-2«) 



(6.21) 2^-^"7r , .,. ^ • 

Now from the definition of (12. 6p . we obtain the right hand side of fl6.13p . as desired. D 

We begin our proof of f l3.14p by applying Lemma \6A] to fl5.40p . with u = s + a, obtaining 

(6.22) M^{k = a,a) = l- Y, ^^. I ^^9.,p,,{s)F{^ + f - «)!,+„ 



a<Y 



2{lay^-Ul -2a- ^s)Res^^,^J-^^^^^^^ds. 

One could wonder why it was beneficial to apply Lemma 16. 1[ The answer is that we wish to 
combine this term with Mr{—1, 1, 1) which has a similar weight function to that appearing 
in (I6.22p . This idea was used in JY]. 

Recall that Cs = \ + e. It is convenient to represent the residue of Ji(2s + 2a, w) at 
w = \ — ahj the value of Ji(2s + 2a, w)/C,{\ + w + a) aX w = \ — a. Then we obtain 

(6.23) M^{k = U,a) = \ J^ ^i / ^^.,/3,7(^W + f " «) 

(a,2l)=l -'^"^f 

a<Y 

Ts+aa'^''^^'Dj^{k = D, a; s)ds. 

where 

Ji{2s + 2a,w) 



(6.24) D^{k = D,a-s)=2F^%,{l-2a-2s)- 

C[2s + 2a)C{^ + w + a) 

Before performing further analysis of this function, we recall from (14. 8 P that 



W=T:—a 



{a,2l)=l " '" -^(^^ 



f-«) 



(6.25) M«(-l,l,l) = i Y: #^/ ^9.,,.is)Fil 

2 ^ ^ ^ a^ 2m J(^) s 

T^+sa^'^+^'Dni-l, 1,1; s)ds. 



a>Y 

2a+2s 



where 
(6.26) 

1 -sp 1 v^ /i(ri)/i( r2)/i(r3) A_a-s,p+s,'r+s{lrir2r3,) 



DJ-lAA;s) 



'R 



^^(2) ,,tt, h-''—'' ^^^^f-,^^ .H^+vf'^^v^^^^ ^W^^^^* 



By analogy with (16. 3p and similar identities in ^, it may not be surprising that the 
following miracle occurs: 

Lemma 6.2. We have 

(6.27) D^i^k = D, a; s) = Dr{-1, 1, 1; s). 
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As in Section I6.H we shall verify this with a computer calculation in Section 16.31 Taking 
it for granted for now, we then have 

(6.28) MM{k = U,a) + MR{-l,l,l) = \ V ^^ [ ^g.,p,,{s)F{l + I - a) 

ia,2i)=i "" ^""'J^^) ' 

Grouping ah into a variable, we see from the Mobius formula that only ab = 1 survives, and 
hence ri = r2 = r^ = 1. Thus 

(6.29) M^(A; = n,a) + M^(-l,l,l) = 

^^^^^— J^^^ -^g.,asW + f - «)r.+.A_,,+,,+.(/)d.. 

By Lemma l4m and Remark [221 we can move the contour of integration to — 1/2 + e, crossing 
a pole at s = only. The residue at s = gives the main term in (I3.14p . and the error term 
is of size 0{X^^'^~^^l^), as desired. 

6.3. Proof of Lemma 16.21 The basic plan is similar to that used in Section [6111 check the 
identity at each Euler factor using a computer. There is a minor issue with convergence that 
we delay explaining for a moment. First consider the case p = 2 (which can be easily done 
by hand). The Euler factor for Dn at p = 2 is 1, while the Euler factor for Dn is 

(6.30) 2(1 - 2-2^-2")(l - 2-i)(l - 2-2^-2")-! = 1. 

Now we add the following definitions to those used in Section [6m w = p2o+2s^ 

(6.31) V{x, y,z) = f: ^""--^n^"'^^'''"^ = E '-^^^^^x'^'^y'z^. 

n=0 ^ a,b,c>0 



and 

n=0 t' a,b,c>0 



(6.32) W{x, y,z) = Y: ''—-^-^;^^ ^ = J^ 1^ 



x^w^yh". 



We view x,y,z,p as free variables; note w = l/{px'^). Then the Euler factor at p\a for D^^ 
is (1 - p-^){l - ^)-\ while for Dr it is 

(6.33) (1 - p"2)(l - wp^'r'll + (1 + P"'r'{-1 + Q'}], 
where 

(6.34) Q' = V -{x + y + z)W + {xy + xz + yz)V - xyzW. 

A computer indeed verifies these are equal. For p \ 2a/, we have that the Euler factor at p 
for Dn is 

(6.35) {l--r\l--){l--){A + Bl 

p w p 

where 

(6.36) A = {1 )-^(l + (l-p-^)[-l + T]), and B = wxU. 
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The corresponding Euler factor at p for Dji is 

ID ID 

(6.37) (1 - p~^){A^^^s,/3+s,^+s{l) + -(1 - -y^[A^a-s,p+s,^+s{l] 

(x + y + z + xyz) ^ ^^^ ^ ^_^ ^ ^_ ^ ^_^ ^ ^ ^ 



-A_a-s,f5+s,-y+sip) + {xy + XZ + yz)A_a-s,f5+s,-y+s{P )]}• 



Vp 

Then observe 

(6.38) ^_„_,,^+,,^+,(l) = 1 + (1 + p-')-'(-l + V), 

(6.39) ^A_,_,,^+,,^+,(p) = (1 + p'Y'W, A_^_s,f,+s,^+s{p') = (1 + p-Y'V. 

VP 

Again, a computer verifies the equahty of Euler factors. We omit the final case p\l (two 
cases, actually, dependong on if /p = 1 or Ip = 2) which follow similar lines. 

Finally, we argue why it suffices to check that the Euler products agree at the desired 
point. Towards this end, we have 

Lemma 6.3. There exists 6 > such that for a,/3,7 small enough compared to 5 we have 
that 

(6.40) ^^("'^^ 



C(v)C(i + W + a)C(| + W + (3)al + w; + 7) 

has meromorphic continuation to the region Re{v) > —6, Re{w) > | — 5. More precisely, 
(I6.40p equals, with (ai,a2,a3) denoting (a,/?, 7), 

.g^^^ Ui<i<j<3 C(^' + 2w + ai + aj) 

C{l + w + v + a)C{l + w + v + ^)C{l + w + v + -f) 

times an Euler product that is absolutely convergent in the above-stated domain. 

Proof. It suffices to consider the behavior of the Euler factors with p \ 2al in which case 

recall the Euler factor is given by (15.241) . This takes the shape 

(6.42) 

(1 -P"^)"Ml + (1 -P~')^^fM^ + 0{p-"'"^n + (1 -p-'')P^^-'"(aa,,,,(p) +0(p"^"2-+^))] 

From this description we can read off f l6.4ip . D 

By Lemmas 16.31 and 14. ![ both sides of (16.271) can be continued meromorphically to a 
domain Re(s) > —6 with some 6 > 0. Thus it suffices to check for each prime p that the 
Euler factor at p for each side agrees with the other side, as claimed. 
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